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Abstract. 

A complete model of the dynamics of scrape-off layer filaments will be rather complex, 
including temperature evolution, three dimensional geometry and finite Larmor radius effects. 
However, the basic mechanism of E x B advection due to electrostatic potential driven by 
the diamagnetic current can be captured in a much simpler model; a complete understanding 
of the physics in the simpler model will then aid interpretation of more complex simulations, 
by allowing the new effects to be disentangled. Here we consider such a simple model, which 
assumes cold ions and isothermal electrons and is reduced to two dimensions. We derive the 
scaling with width and amplitude of the velocity of isolated scrape-off layer filaments, allowing 
for arbitrary elliptical cross-sections, where previously only circular cross-sections have been 
considered analytically. We also put the scaling with amplitude in a new and more satisfactory 
form. The analytical results are extensively validated with two dimensional simulations and 
also compared, with reasonable agreement, to three dimensional simulations having minimal 
variation parallel to the magnetic field. 


Filaments are a prominent feature of the scrape-off layer (SOL) of tokamaks, in L-mode 
and in H-mode both during and between ELMs[l], and in other magnetised plasmas. They 
provide a significant component of the particle transport [2], especially in the far SOL, and 
so may have a strong impact on particle fluxes to the first wall and divertor. Data from gas- 
puff imaging diagnostics [3] suggests that filaments may sometimes have significantly elliptical 
(rather than circular) cross-sections. Most of the theoretical work on SOL filaments has been 
done in simplified, two dimensional models, see [4] and [5] for reviews, but recently there has 
been an increasing amount of attention given to more realistic models, for example using three 
dimensional simulations[6-8] or including finite Larmor radius effects[9]. These developments 
have motivated us to reexamine and extend the earlier analytical work to try to give as complete 
a physical picture of the two dimensional mechanisms that regulate filament motion as possible, 
in order to facilitate the interpretation of more complicated models applied to isolated filaments, 
and ultimately to SOL turbulence siniulations[10-12]. 

Our subject here is the scaling of filament velocity with various parameters of the filament, 
in the two dimensional limit where parallel variation may be neglected and assuming cold ions 
and isothermal electrons. This represents the very simplest model that can capture the basic 
mechanism of filament motion. By characterising quantitatively the basic physical processes 
driving filament motion in this simple model, we will be able, when analysing more complicated 
models, to identify the deviations from this behaviour. These deviations can then be ascribed 
to the extra physics in, for example, three dimensional or non-isothermal simulations. The aim 
therefore is to provide a theoretical tool to aid the interpretation of more complicated models. 


as part of a programme of research building systematically towards models which are both 
well enough understood to be trusted and also realistic enough to be quantitatively compared 
to experiment and used predictively for future machines. The particular model we consider 
contains two mechanisms that limit the filament velocity: inertia and sheath currents. While 
inertia is universally present, the sheath current (as modelled here) is only relevant on the 
assumption that parallel resistivity is small so that parallel currents may reach the sheath 
unimpeded. Where this assumption is invalid, whether due to cold plasma[13] or interaction 
with neutrals in the divertor or to large magnetic shear near X-points, different mechanisms will 
come into play; for some examples see [4]. If it is desired to include other effects, as different 
closures for the two dimensional equations, we hope that the framework of the calculation 
provided here will make further additions relatively straightforward. 

We present here a calculation of the scaling of the maximum filament velocity allowing 
arbitrary elliptical cross-sections, which have not previously been considered analytically. We 
also aim to clarify the physical mechanisms while giving more detail than has previously been 
published. In particular, the scaling with amplitude, though considered in [14] and mentioned 
briefly in several works[15-17], has not previously been given a satisfactory analytical treatment; 
however, [18] has provided a complete description of the amplitude scaling in the inertial regime 
from simulations, highlighting the importance of not making the Boussinesq approximation. We 
use two dimensional simulations of isolated filaments both to identify the qualitative features 
that go into the calculation, which is especially important for the inertial regime as it is strongly 
non-linear, and also to validate comprehensively the scaling laws obtained. 

1. Velocity scaling 

1.1. Physical picture 

An isolated filament is, under the assumptions here, a density fluctuation comparable to 
or larger than the background in amplitude, extended along the magnetic field and with a 
monopolar structure perpendicular to it. The basic mechanism of filament motion is that there 
is a current source due to the inhomogeneity of the magnetic field which propels the filament, 
balanced in the simple model considered here by inertia of the surrounding background plasma 
and by dissipation from currents through the sheath where the magnetic field intersects the 
wall at the targets. This description is in some respects inspired by the ‘equivalent circuit’ 
picture[4, 19]. In contrast to the qualitative description in [4], here we consider drifts/currents 
only in the fluid model, rather than considering the particle drifts. In order to obtain simple 
expressions for the filament velocity, we reduce the problem to two dimensions by assuming zero 
variation parallel to the magnetic field and look only for a maximum velocity of the filament. 
First we sketch the physical picture, before turning to the full calculation in Section 1.2. 

Drive Due to the curvature of the magnetic field, represented by k = —x/Rc where Rc is 
the radius of curvature, and the change in its magnitude with radius, VB = —Bx/Rc, the 
diamagnetic current Jjia = B~^V{nT) x b {b = B/B being the unit vector in the direction 
of the magnetic field) has a non-zero divergence, as sketched in Figure 1. This divergence 
gives a current source whose strength is proportional to the pressure gradient in the direction 
perpendicular to k, and VB, which we will call z. In order to maintain quasineutrality, this 
current must be closed. Depending on the parameters of the filament, this may occur either 


through the polarisation current (leading to inertial evolution) or the parallel current (giving 
rise to sheath current dissipation). 
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Figure 1. Diamagnetic current in a filament: The width of the line represents the magnitude 
of the current (integrated in the 6-direction). The diamagnetic current (left) can be split into 
a divergence-free part (centre) and the remainder (right). This last, divergent part must be 
closed by parallel or polarisation currents and so provides the drive for filament motion 


Inertial evolution In the absence of dissipation mechanisms (viscosity or sheath currents), 


there can be no steady state, since the energy of the system increases monotonically[18]. 
Nevertheless, the filament reaches a maximum velocity, even without dissipation. This contrasts 
with hydrodynamic drag on a solid object, where the velocity is ultimately limited by viscous 
dissipation either in a surface layer around the object or in the wake behind it. The mechanism 
causing the hlament velocity to saturate is that it expands as it moves, drawing in fluid from 
the ‘background’ flow. This is rather like a rocket engine in reverse; instead of generating thrust 
by expelling material, the filament resists acceleration by absorbing material. The maximum 
velocity occurs when the rate at which momentum is acquired through this absorption balances 
the J X B force from the polarisation current. 

Sheath current dissipation When the current source closes through the parallel current it is 
regulated by the sheath boundary condition, (12), which determines the current that passes 
through the Debye sheath to the target plates. For small currents the sheath behaves as an ideal 
resistor, with current proportional to the electrostatic potential, 0. In steady state the potential 
is then proportional to the current source, behaving as a simple electrical circuit (Figure 2). 
The motion of the filament is then just the E x B drift due to this potential. 



Figure 2. Equivalent circuit for a filament in which the diamagnetic current closes through 
the parallel current path 
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Figure 3. Density (colour-map) and potential (contours; solid for (j) > i), dashed for </) < 0) for 
a filament in the sheath regime, at maximum T4 


1.2. Scaling calculation 

We use a coordinate system that takes x in the ‘radial’ direction (anti-parallel to the magnetic 
field curvature and VB), y along the magnetic field and z in the binormal direction. We take 
the filament to be a monopolar density structure (possibly distorted, in the inertial regime), 
characterised by two length scales, 6^ and 6^ in the x- and z-directions, and an induced, dipolar 
potential structure. To aid visualisation of the configuration, the density and potential from a 
simulation (see Section 2) of a filament in the sheath current regime are shown in Figure 3. We 
assume cold ions and isothermal electrons, neglect electron inertia and use Bohm normalisation 
(as in [8], times are normalised to the ion cyclotron frequency flci = lengths to the 

hybrid Larmor radius ps = Cg/ffci, where Cg = is the sound speed, and electrostatic 

potential by Tg/e, with e the elementary charge, Bq the magnetic field strength, m\ the ion 
mass and Tg the electron temperature; in addition densities are normalised to a typical value, 
usually the background density). 

The system is driven by the diamagnetic current, and its behaviour determined by whether 
the divergence of the diamagnetic current closes principally through the polarisation or the 
parallel current. We now consider each current in turn. 


Diamagnetic current For cold ions, the diamagnetic current is due just to the electron 
diamagnetic velocity. Neglecting electron inertia and viscosity, the electron momentum equation 

is 

V.xB 


0 = —Vn -|- nV0 -|- n- 




( 1 ) 


Taking the cross product with h gives the perpendicular electron velocity with two contributions, 
the E X B and diamagnetic velocities: Ve = Vexb + Vdia, with 
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Here we must not yet take B = Bq (Hq being the constant magnetic field used for normalisation) 
until we have taken its gradient in the divergence of the diamagnetic current: 
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Figure 4. Density (colour-map) and streamlines of flow velocity in the frame of the filament 
(white) at the time when maximum velocity is reached 
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Figure 5. Vorticity (colour-map) and streamlines of flow velocity in the frame of the filament 
(white) at the time when maximum velocity is reached 
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for a toroidal magnetic field B = e^/R at a major radius R = Rq [Rq denoting the radius of 
curvature of the magnetic field and being the unit vector in the toroidal direction) much 
larger than the size of a filament and defining g = gx = x/Rq. 

Polarisation current Define the filament velocity, V{, as the velocity of the centre of mass of 
the perturbation, 

V{ = J dPx{n-no)V ( 4 ) 

and the ‘filament frame’ as that moving at with respect to the lab frame. In the filament 
frame there are two distinct regions (see hgure df): the ‘exterior region’ outside the filament 
where the velocity is always to the left (the streamlines are open); and the ‘filament region’ 
where the velocity circulates (the streamlines are closed). 

f The figures 4, 5, 6 and 7 in this section are taken from a simulation with filament amplitude eight times 
the background using the equations (17), (18) and (19) but, in contrast to Section 2, the curvature terms 
(proportional to g) in (17) are dropped and sheath currents and viscosity are dropped from (18) (i.e. Ly —^ oo 
and /Ti = 0). This ensures that we are examining here a purely inertial filament without E x B compressibility, 
which we do in order to understand the qualitative features of this regime. 







Figure 6. Evolution of three density contours with time, hz is the height of the contour: the 
difference between the maximum and minimum of z on the contour. The snapshots on the 
right show the contours (at levels n = no (l + A/e*) for i = 1,4,9 from inside to outside, 
corresponding to the lines in the left plot) at the times indicated by the dashed vertical lines, 
t = 87 is the time at which the filament reaches maximum velocity. The horizontal dotted lines 
indicate the upper and lower extents of the contours at t = 0 


There is always a force on the hlament from the polarisation current Jpoi x b force, which 
tends to increase its momentum (in the x-direction). For the velocity to be constant this gain in 
momentum must be balanced. In the frame of the background plasma, the balance is provided 
by absorbing stationary background plasma, which increases the mass of the hlament (without 
adding extra momentum), allowing the velocity to remain constant. Alternatively, considered 
in the frame of the hlament, the background plasma has negative x-momentum relative to the 
hlament, providing a loss of momentum to balance the gain due to the Jpoi x b force. 

We need to evaluate the rate at which the hlament gains momentum by absorbing plasma 
from the background. E x B how is (approximately) incompressible, so areas that are advected 
in the how, in particular the areas of density contours, are preserved. Therefore an increase in 
the area of the hlament can only be due to its absorbing plasma from the background (due to 
changes in the how pattern so that the closed streamline region, the ‘hlament region,’ grows). 

Even hlaments initialised with a Gaussian density prohle (as in the simulations described 
below) generate during their evolution a much sharper boundary between ‘hlament’ and 
background than was present in the initial (Gaussian) density prohle, by the time they reach 
maximum velocity. As we see in Figure 6, the outlying parts of the initial density distribution 
are either gathered into the hlament, or are left behind altogether (although only a negligibly 
small fraction of the total density perturbation is left behind). The result is that by the time the 
hlament reaches its maximum velocity the heights of the contours of a wide range of densities 
(several orders of magnitude) are gathered together close to a single value, just inside the closed 
streamline marking the boundary of the hlament (see Figure 4), and the heights are increasing 
at the same rate. Therefore the plasma absorbed from just outside the hlament boundary has 
a density (almost) identical to the background density, hq. 

We denote the stream function in the hlament frame by (whose contours are the 
streamlines plotted in Figures 4, 5 and 7). V’ is related to the electrostatic potential (which 
is the stream function in the lab frame) by 


'0 = 0 — ViZ. 


( 5 ) 













Figure 7. The stationary points of the flow (the dotted lines are streamlines of flow velocity 
in the frame of the filament) pin the ends of density contours which are ‘inflated’ by the flow 
between the stationary points, with velocity ~ Vc 


In the exterior region the density gradients and vorticity are negligible, as we see in figures 4 
and 5, which means that = 0 and hence = 0- The boundary conditions are constant 
-0 on the filament boundary (since we define the boundary by a streamline) and V = —V{ 
asymptotically at large distances from the filament. The flow in the exterior region is thus the 
solution to a linear problem and scales linearly with the magnitude Vf. The average speed of 
the plasma absorbed by the filament from the external flow therefore scales like Vf. 

Putting these two results together, the plasma captured as the filament increases its area 
will have an average x-momentum density 


(Px) ~ -noVf 


up to some geometrical factor depending on the shape of the filament, which we assume does not 
change much, the assumption being justified a posteriori by the fact that the result obtained 
describes the simulations well. 

The rate at which the filament gains momentum from the background is then 


dP. 

dt 


dA 


-noVfL|| 




where A is the area of the filament (the area enclosed by the largest closed streamline). Finally 
we need to estimate dA/dt. The filament increases its area because the crescent-shaped density 
contours are ‘inflated’ (since they cannot be compressed, their areas being preserved) by the 
flow on the axis, between the stagnation points of the circulating flow within the filament, 
which are separated by a distance ~ 6z, as illustrated in Figure 7. If the average velocity (in 
the filament frame) of that flow is Vc then 




Vc6z 


( 6 ) 


The only velocity scale in the system is Vi, so (in the absence of any more exact calculation), 
we estimate K Vf so that 


dP. 

dt 




(7) 


The total force on the filament is in the positive x-direction, with magnitude {P denotes 
the filament region and dP its boundary) 
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integrating by parts twice. We may neglect the surface terms because firstly d^x Jpoi,z ~ 0 
in the exterior region, since the flow pattern is almost constant and so the external plasma 
must be close to global force balance. Since also V ■ Jpoi = 0 in the exterior region 
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and Jpoi vanishes at inhnity (r^Jpoi —)■ 0 as r —>■ oo since otherwise the entire background 
plasma would accelerate). Secondly (n — uq) Ri 0 on the boundary, so f dx [zg {n — ~ 0, 

where z^ and denote the upper and lower positions of the boundary at a particular x. The 
area integral is just the total number of particles in the perturbation, which is conserved from 
the initial state (and almost all contained within J”) and which we may estimate as 


IF 


d‘^x {n — no) ~ AuQSxdz 


( 10 ) 


where A is the amplitude of the hlament relative to the background and A ~ d^dz is the area of 
the hlament (we assume that the changes in A with time are small compared to its magnitude). 

Therefore the maximum hlament velocity is attained when the rate of momentum gain 
from the exterior plasma, (7), balances the total force, (8), so 

Vf ~ \/gAdx 


Parallel current Linearising the sheath boundary condition. 


^11 
" sheath 

and integrating over y, 


= insheath (l - e ~ ±n(j) 


dyV ■ J\\ ^ 2n(f). 


(11) 


( 12 ) 


(13) 


Since the hlament velocity is approximately the central E x B velocity of the hlament (in the 
lab frame), we can further estimate V{ ^ (j)/dz and hence 


dyV ■ J\\ ndzVf. 


(14) 


So balancing this with the diamagnetic divergence, (3) (integrated over y), the velocity in this 
regime is 

The interpretation of the factor (n — no) /n ~ Aj (1 + /3A) and the value of the constant fJ will 
be discussed in Section (2.2). 






We have two regimes for filament dynamics, the inertial regime with maximum velocity 
given by (11) for narrow filaments and the sheath current regime with velocity given by (15) 
for wide filaments. Interestingly, the velocity depends on 5^ but not 5z in the inertial regime 
and conversely on 6z but not 5x in the sheath current regime. The latter point was noted in the 
seminal paper [20], in the case without background plasma where a separable solution exists, 
but seems to have been neglected since then. The transition between the two regimens occurs 
when the velocities (11) and (15) are comparable, so 

( 16 ) 


2. Comparison with simulations 


We have used two dimensional simulations to validate the scalings derived in Section 1.2. The 
equations used, based on those in [8], represent a filament assumed to have negligible variation 
along the magnetic field; closure is given by integrating the three dimensional system over the 
parallel direction to give equations for the density and vorticity 
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(18) 


11 = V ■ {nV±<l>) 


(19) 


with d/dt = {d/dt + b ■ V0 x V^, V exb = bx V0 and, as before, g = gx. (18) may be derived 
by taking the parallel component of the curl of the ion momentum equation 

= Jpoi X h + /iiVi (nVO (20) 

and replacing V ■ Jpoi using quasineutrality, V ■ Jdia + V ■ Jpoi + V ■ Jy =0. (19) is inverted 
for 0 using a new multigrid solver currently under development for BOUT++, to be reported 
elsewhere. The dimensionless parameters used are L\\ = 11000, g = 2.5 x 10“^, fin = 1-5 x 10“® 
and Pi = 4 X 10“^; for this theoretical study the dissipative parameters have been reduced by 
a factor of 100 compared to those used in [8], based on the expressions in [21], in order to 
minimise their effect on the inertial regime, since we do not treat viscosity here. Filaments are 
initialised as Gaussian density fluctuations with elliptical contours, tilted at an angle a to the 
;r-direction, on a constant background, 

n(t = 0) = no ^1 + Aexp (21) 

where A is the amplitude of the filament, e is the ratio of the lengths of the axes of the ellipse, 6 
is the geometric mean of the lengths of the axes, x' = x cos a-hz sin a and z! = z cos a — x sin a. 
The diagram in figure 8 shows the configuration with the corresponding length scales, in the 
X- and z-directions respectively. 


dx = d 
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( 22 ) 




( 23 ) 
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Figure 8. Schematic showing sizes and orientation of a tilted, elliptical filament cross-section 




XXX 

6=5.0, a 

= 90.0° 

)->->- 6=3.0, a 

= 60.0° 

1 1 1 6=1.0, a 

o 

O 

o 

II 


6=3.0, a 

o 

d 

II 

CCO e=5.0, a 

II 

o 

d 

o 

+++ 

6=5.0, a 

= 60.0° 

-<-<-( 6=2.0, a 

= 90.0° 

o o o 6=2.0, a 

o 

q 

o 

II 

□□□ 

6=3.0, a 

= 30.0° 

OOo 6=5.0, a 

= 30.0° 

YYy 

6=3.0, a 

= 90.0° 

— 6=2.0, a 

= 60.0° 

7W 6=2.0, a 

= 30.0° 







Figure 9. Inertial scaling vs. Sx (left) and sheath cnrrent scaling vs. Sz (right) with A = 4 for 
several values of e and a 

which appear in the velocity scaling above. Filament velocities are measured as the maximum 
velocity in the ^-direction of the centre of mass of the density above background, i.e. n — no- 
These simulations, as well as the three dimensional ones in Section 3, have been implemented 
using BOUT++[22, 23]. 

2.1. Scaling with size and shape 

The results of a scan in width S, ellipticity e and inclination a, at constant amplitude A = 4, 
are shown in Figure 9. For four different e and four different a, scanning 6 results in a region 
of scaling, corresponding to the inertial regime, and a region of scaling corresponding 
to the sheath current regime. Moreover, although we have a three parameter space, {5, e, a}, of 
filament sizes and shapes, a single combination (h^, in the inertial regime and Sz in the sheath 
current regime) entirely determines the filament velocity, at a particular amplitude (the effect 
of amplitude will be discussed in Section 2.2); this is evident from Figure 9 since all the points, 
for any e and a (in the region where they follow y/6^ or 5“^ scaling) have not only the same 
gradient but also the same absolute magnitude. 

The results for e = 3, a = 90° and e = 5, a = 90° exceed the scaling for Sz ~ 20 — 50. 
This is due to the non-linearity of the sheath boundary condition which is included in the 
simulations, but not in the scaling. The maximum value of the potential is 0 ~ 1.0 at Sz ~ 26 
for e = 3 and (j) ~ 1.0 at Sz ~ 34 for e = 5, so it makes sense that the deviation of exp (—0) from 
its linearised form is noticeable here. The effect of the non-linearity is to decrease the magnitude 
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Figure 10. Amplitude scaling in the inertial regime, 5^ = 5^ — 0.5, (left) and in the sheath 
current regime, Sx = Sz = 200, (right). The solid red line shows the fitted ansatz, (24), and the 
dashed black line shows best fit line oc -x/A 


of the negative lobe of potential and increase the positive lobe. However, the enhancement of 
the positive lobe must be larger in order to allow the same magnitude, jo, of current through 
the sheath: |0_|_| / |0_| = — In (1 — Jq) / In (1 + jo) > 1. The overall effect is therefore to increase 
the hlament velocity, since the increase from the positive lobe outweighs the decrease from the 
negative lobe, resulting in a larger Ez- 

The trend lines are V ~ 5.7 x for the inertial scaling and V ~ 11^7^ for the sheath 

scaling. These intersect at ~ 8.2 whereas (16) suggests a crossover at ~ 11.9 

for these parameters (using fj from Section (2.2)), in good agreement up to the order unity 
factors which are not hxed by analytical scaling arguments. 


2.2. Amplitude scaling 


The amplitude scaling predicted for the inertial regime is simple, just V oc y/A, where the 
maximum fluctuation amplitude is H = max [(n — no) /no\. For the sheath current regime the 
prediction is a little more complicated. We may estimate Vz{n — uq) ~ Anojbz, but we need 
to hnd a representative value of n to estimate z{n — no) in (15). We hnd, by htting the 
simulation results, that the appropriate value is neither the maximum, n ~ (1 + A)no nor the 
background n ~ no, but rather an intermediate value n ~ (1 + (jA)no, where 0 < /3 < 1. /I 
corresponds to some point on the density prohle at a distance 5z'. (3 = exp(—Then Xj as 
a function of A with all other parameters held constant is 


Hf oc 


A 

1 + (3 A 


(24) 


if the representative value of n is found at the same relative position, 6z/Sz for any A. This does 
indeed seem to be the case, as illustrated in Figure 10. For an amplitude scan with = Sz = 0.5, 
which is in the inertial regime V oc y/A hts well, in agreement with [18], with e ~ 0.016, where 
£ is the root-mean-square relative error. For an amplitude scan with S^ = Sz = 200, which is in 
the sheath current regime, (24) is a good £t, e ~ 0.014, and gives f3 ~ 0.31. These conclusions 
also hold when varying the ellipticity, see Appendix A. 

Comparing this to previous work: Angus and Krasheninnikov[18] established numerically 
the y/A scaling in the inertial regime, which we have tried to explain analytically above; Kube 
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Figure 11. Comparison of two dimensional (lines) and three dimensional (markers) results for 
A — 2 and several e 


and Garcia[14] consider the amplitude dependence in some detail, where the scaling analysis 
starts from the vorticity equation and finds similar results to those given here, except that 
they use the Boussinesq approximation in their simulations and so do not hnd -s/A in the 
inertial regime, while in the sheath current regime they assume a form equivalent to setting 
= 1 which prevents them from finding a quantitative, analytical amplitude scaling (their 
fit coefficients vary with amplitude, whereas here we have only a single, constant coefficient, 
/3, to be derived from simulations); Angus et ah [6] give, albeit briefly, a derivation starting 
from the momentum equation, but neglect the amplitude dependence of the drive; Theiler et 
ah [17] in contrast use an interchange instability growth rate to estimate d/dt ~ 7interchange, and 
hnd a linear scaling with amplitude of the hlament velocity in the inertial regime (which is 
contradicted here, emphasising that it is the non-linear, advective, term that is relevant); their 
derivation in [17] follows Garcia et ah [15, 24], but there the ‘ideal interchange rate’ is dehned, 
without explanation, as 


7 = 


(Iff 


(where g corresponds to our g, I to Ly, A6' is the huctuation amplitude and 0 is the background 
density) which includes an amplitude dependence giving a square-root scaling of the hlament 
velocity in the inertial regime. 


3. Three Dimensional validation 

For reasons of computational expense, the three dimensional simulations in this section were 
done using the Boussinesq approximation (as were the two dimensional simulations just in this 
section, for consistency). This approximation is known to give the wrong amplitude scahng[18], 
but seems not to affect much the scaling with size and shape (at least the -y/h]) inertial regime 
and sheath current regime may be found either with or without the approximation). Here 
the normalised background density in the two dimensional simulations was set to uq ~ 1.48 to 
be consistent with the (source-driven) background used for the three dimensional simulations, 
which normalise to the equilibrium density at the sheath entrance. The amplitude was set 
to a lower value than in the two dimensional case, A = 2, to avoid the possibility of drift 







wave instabilities playing a role in the three dimensional simulations. For more details on the 
equations used in and implementation of the three dimensional simulations, see [8]. The three 
dimensional filaments were initialised with no variation in the parallel direction, in order to 
correspond to the two dimensional calculation as closely as possible; investigation of really 
‘three dimensional’ effects due to parallel gradients is left for future work. Figure 11 shows that 
the three dimensional simulations follow the two dimensional trends closely. 

The velocity of three dimensional filaments in the sheath current regime is slightly larger 
than the corresponding two dimensional ones. We attribute this to the variation in the 
background density: since this decreases near the sheath, the potential needed to drive the 
same sheath current is slightly increased and it follows that the filament velocity must then 
increase slightly. 

The bump on the e = 1/3 curve between 6 = 20 and 6 = 30, the part most affected by the 
non-linearity of the boundary condition (see Section 2.1), is due to the eventual fragmentation 
and subsequent acceleration of the filaments. In other words, it is due to their ceasing to travel 
as coherent structures and so its investigation is a subject beyond the scope of this paper. 

4. Discussion 

A number of works[16, 17] have analysed filament motion by analogy with the theory of linear 
instabilities; this has been termed the ‘blob correspondence principle’ [16, 25]. Our analysis above 
shows that the picture is not quite so simple; the nature of the filaments as coherent, non-linear 
objects is important. Even in the sheath current regime, the density and potential fields have 
qualitatively different structures, being respectively monopolar and dipolar (Figure 3), rather 
than the identical structures, up to a phase shift, that they would have in linear theory; in 
the inertial regime the difference in the form of the fields from a simple mode structure is 
even more marked (Figures 4 and 5). Thus the nature of the physical processes involved in 
limiting hlament velocity is clearer when considering the filament itself instead of analysing the 
governing equations by analogy or ‘correspondence’ to linear theory. 

It is notable that in the inertial regime rather fine structures are formed, especially in the 
vorticity (see Figure 5). This suggests that both finite Larmor radius effects and viscosity may 
have significant effects on filaments in the inertial regime, which could provide an interesting 
topic for future study, perhaps using a gyrofluid model. 

5. Conclusions 

We have given here the first calculation of SOL filament velocity to include the effect of 
the filament shape and have also clarified the role of the filament amplitude. The analytical 
scaling calculations have been extensively validated by two dimensional simulation results, and 
also compared (with good agreement) to three dimensional simulations in which the filaments 
are initialised without parallel variation. Thus we now have a complete understanding of the 
mechanisms of filament propagation in the simple limit considered here. This understanding 
provides a solid foundation for the interpretation of filament motion in more complicated, more 
realistic models. 
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Appendix A. Amplitude scaling fits 

In order to verify that the form of the amplitude scaling is independent of the shape of the 
hlament, we repeated the scans described in Section 2.2 for several different ellipticities, keeping 
6x constant for the hlaments in the inertial regime and 6z constant for the hlaments in the sheath 
current regime so that (according to (11) and (15)) the velocities of the hlaments should be 
independent of the ellipticity. 

The root-mean-square relative error, e, is dehned by = A ; where V) is the 

velocity measured from the simulation with amplitude Ai, V{A) is the predicted scaling and N 
is the number of simulations in the scan. 


<5. 

£ 

0.5/5 

0.052473 

0.5/3 

0.039863 

0.5/2 

0.029693 

0.5 

0.015660 

0.5 X 2 

0.022420 

0.5 X 3 

0.031436 

0.5 X 5 

0.044025 


Table Al. Inertial regime: Comparison of e for the scaling V oc \fA at = 0.5 for several 
values of Amplitudes used for the analysis were A = 1,1.5, 2,3,4, 5, 6, 7,8. 












dx 

P 

£ 

200/5 

0.30985 

0.014128 

200/3 

0.30994 

0.014123 

200/2 

0.30997 

0.014123 

200 

0.30999 

0.014121 

200 X 2 

0.30999 

0.014120 

200 X 3 

0.30999 

0.014121 

200 X 5 

0.31000 

0.014122 


Table A2. Sheath current regime: Comparison of (i) values of f3 inferred by least squares 
regression on the relative error between the ansatz (24) and filament velocities measured 
from simulations and (ii) values of £ for this regression, at 5^ = 200 for several values of 
5x- Amplitudes used for the analysis were A = 1, 1.5, 2,3,4, 5, 6, 7,8. 
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